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Abstract 

We consider the interaction between the Hermitian world, repre- 
sented by a real delta- function potential — aS(x), and the non-Hermitian 
world, represented by a PT-symmetric pair of delta functions with 
imaginary coefficients i(3(5(x — V) — 5{x + L)). In the context of stan- 
dard quantum mechanics, the effect of the introduction of the imagi- 
nary delta functions on the bound-state energy of the real delta func- 
tion and its associated wave-function is small for L large. However, 
scattering from the combined potentials does not conserve probability 
as conventionally defined. Both these problems can be studied instead 
in the context of quasi-Hermiticity, whereby quantum mechanics is 
endowed with a new metric rj, and consequently a new wave-function 
^(x), defined in terms of the original wave-function ip(x) by means of 
ij. In this picture, working perturbatively in (5, the bound-state wave- 
function is actually unchanged from its unperturbed form for \x\ <C L. 
However, the scattering wave- function, for |x| S> L, is changed in a sig- 
nificant manner. In particular there are incoming and outgoing waves 
on both sides of the potential. One can then no longer talk in terms 
of reflection and transmission coefficients, but the total right-moving 
flux is now conserved. 
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1 Introduction 



Since the resurgence of interest in non-Hermitian Hamiltonians initiated by 
the paper of Bender and BoettcherpQ the subject has gone through several 
stages. First there was the search for soluble non-Hermitian Hamiltonians 
with real spectra (see, e.g. [2]). The next stage concerned the introduc- 
tion of a new Hilbert-space metric in order to obtain positive probabili- 
ties and so regain a physical interpretation of the theory. This was first 
done in PT-symmetric theories [3] by introducing the so-called C operator 
to form the metric CPT. A more general framework was formulated by 
Mostafazadehjl], which among other things established the connection to 
earlier work by Scholtz et al. [5] and showed[H] that such a Hamiltonian H 
was related by a similarity transformation to an equivalent Hermitian Hamil- 
tonian h. Subsequent work showed how the metric rj could be constructed, 
sometimes exactly, but more typically in perturbation theory for a variety of 
models [3 El i- 

More recently some attention has been given to situations where a non- 
Hermitian system interacts with the world of Hermitian quantum mechanics. 
For example, Ref. [10J examined a non-Hermitian analogue of the Stern- 
Gerlach experiment in which the role of the intermediate inhomogeneous 
magnetic field flipping the spin is taken over by an apparatus described by 
a non-Hermitian Hamiltonian. This type of set-up has been further elab- 
orated by Assis and Fringjll] and Guenther al.[l2]. and Most af azadeh [13] 
has emphasized that the effect relies on non-unitarity in some guise or an- 
other (see also p3])- The subject has also taken an experimental turn in the 
form of optical lattices whose refractive index can be tailored to make them 
PT-symmetric (see, for example, |15j). 

An earlier attempt at understanding the conceptual issues involved in 
scattering from a non-Hermitian potential was given in Ref. [16]. There we 
primarily considered a single delta-function at the origin with a complex 
coefficient. Here unitarity is not conserved within the framework of conven- 
tional quantum mechanics, but if one instead constructs the metric i] and the 
corresponding transformed wave-function \l/ (see Eq. (|2"U1) below) unitarity is 
restored but at the price of a drastic change in the physical picture, whereby 
there are now incoming and outgoing waves on both sides of the potential. 
A somewhat less drastic, but still significant, change is found by ZnojilfTT] 
in a discretized model of scattering, where the metric does not mix up in- 
coming and outgoing waves, but instead changes the normalization of the 
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flux on either side of the scattering centre. In Ref. [16J we touched on the 
potential model that we address in the present paper, but only in the context 
of conventional quantum mechanics. On the basis of those calculations we 
speculated that the treatment of bound-state problems should be essentially 
unaffected by the introduction of distant non-Hermitian scattering poten- 
tials, and that it is only when the system physically interacts with those 
potentials that a drastic change in the formalism is required, if indeed the 
non-Hermitian potentials are regarded as fundamental rather than effective. 

In the present paper we return to that model, which we are now able to 
treat in the quasi-Hermitian picture as well, thanks to a general prescrip- 
tion due to Mostafazadeh[T8] for calculating r\ perturbatively in the case of 
a superposition of delta functions. The model is introduced in Section 2, 
where it is treated in the framework of conventional quantum mechanics. As 
was already seen in Ref. [16] . the bound-state energy and wave-function are 
affected only by exponentially small terms when L — > oo, but the scatter- 
ing does not preserve unitarity, with R + T ^ 1. In Section 3 we instead 
treat the problem in the quasi-Hermitian framework, calculating the metric 
to first order in (3 but to all orders in a. For the bound state we find that 
the new wave-function \1/ is actually equal to the old wave-function ip in the 
absence of the non-Hermitian part of the Hamiltonian. For the scattering 
problem we find the same type of wave-function that we previously found for 
the single complex delta function, whereby \1/ contains an incoming wave on 
the right-hand side, in contrast to ip. Finally, in Section 4, we discuss the 
significance of these results. 

2 Delta-Function Model 

The Hamiltonian we analyze in this paper has a potential consisting of three 
delta functions: 

H = p 2 - a5(x) + i(3(5(x - L) - 8(x + L)). (1) 

The first component of the potential, — aS(x), is a delta function based on 
the origin with a real, negative coefficient. When (3 = this Hermitian part 
of the potential supports a single bound state with energy E = — k 2 , where 
k = |a. The second component proportional to f3 consists of two delta 
functions based at x = ±L with imaginary coefficients ±i/5, designed to be 
PT-symmetric and have real energy eigenvalues. What is at issue is how the 
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introduction of this non- Hermit ian piece of the Hamiltonian affects both the 
scattering wave-function and the bound-state energy and wave-function. Of 
particular interest is the case when L is large, i.e. when the non-Hermitian 
pieces are distant from the Hermitian potential based at the origin. 



2.1 Bound State 



Let us start by setting (3 = 0. As already stated, there is then a single bound 
state, with (unnormalized) wave- function 



-k\x\ 



where k = \a in order to satisfy the continuity condition 

^'(0 + )-V'(0_)=a^(0). 



(2) 



(3) 



Let us now repeat this calculation for (3 ^ 0. The bound-state wave- 
function then has the (PT-symmetric) form 



We KX 

Ue KX + Ve~ KX 

JJ* e -nx + y* e nx 

W*e~ KX 



x < -L 
-L < x < 
< x < L 
L < x 



(4) 



with^(— x) =ip*(x). Applying the continuity conditions ip{— L + ) = ip{— L_), 
ip'(-L + ) - il/(-LJ) = -i(3il){-L) at x = -L we find 



U 



[1-0)W 



V = i[3e-" KL W, 

where f3 = (3/ (2k) . The continuity conditions at x = are 

U + V = U* + V* 
k[(U + U*)-(V + V*)} = a(U + V). 

The first of these gives 



(5) 



(l-i(3) +i(3e~ 2KL W 
which can be satisfied by taking 

W — (1 + 0) — i{3e~ 2KL 



(1 + 0) - iPe~ 2KL } 



W* 



(6) 



(7) 



(8) 
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up to an overall real normalization constant. On substitution into the second 
continuity condition we find, after some algebra, the eigenvalue equation 



-2kL\2 



2k l+(3 2 (l-e- 4KL ) =a l + (3 2 (l-e 
which, after some additional manipulation, can be written as 

(l + a)/3 2 e- 2Ki (l-e- 2KL ) 



a 



1 + 



1 + (3 2 (1 -e 



-2kL\ 



(9) 



(10) 



where a = a/ (2k). In this form it is clear that (i) k — > |a as L — > oo, and (ii) 
the first correction to k is of order f3 2 . It is interesting that the wave-function 
does not depend explicitly on a, only through the relation between k and a. 



2.2 Scattering Wave-Functions 

Again let us start with (3 = 0. The scattering wave-function, in the situ- 
ation where a plane wave comes in from the left and is either reflected or 
transmitted at the delta-function potential, is of the form 

^)^{^r +Be " ib :t° (id 

The coefficients A and B are determined by the continuity condition ([3]) as 

A = 1 — ia 

B = ia, (12) 
where a = a /(2k), giving reflection (R) and transmission (T) coefficients 
T= 1/\A 2 \ =1/(1 + a 2 ) 

(13) 

R= \B 2 \/\A 2 \ =a 2 /(l + a 2 ). 

Since the potential for (3 = is real, the scattering is unitary, with R + T = 1. 
For (3 the wave-function has the general form 

Ae ikx + Be~ ikx x < -L 

. CtP' + De-** -L<x<0 

' ^ Ee ikx + Fe~ ikx < x < L y } 

e lkx L < x 
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Working from the right, we first apply the continuity conditions at x = L, 
to obtain 

E = 1-/3 

F = pe 2ikL , (15) 

where (3 = (3/ (2k). Then, applying the continuity conditions at x — 0, we 
obtain 

C = (1 - za)(l - $) + ia(3e 2ikL 

D = (l + ia)(3e 2ikL + ia(l- $). (16) 

Finally, applying the continuity conditions at x — —L we obtain, after some 
algebra, 

A = (l-(3 2 )(l-ta)-2ta(3 2 e 2ikL + (3 2 (l + ia)e ilkL (17) 
B = /3(1- + ta)e 2ikL ~ (1 ~ ta)e~ 2ikL ] + t&0 2 + (1 - /3) 2 ]. 

Note that A only involves /3 2 , whereas B involves (3 linearly. The coefficients 
are at most linear in a. As expansions in (3 we have 

A = l-ta + 0{(3 2 ) (18) 
B = ia + 20[sm 2kL -2a sin 2 kL] + 0((3 2 ). 

Thus the transmission and reflection coefficients are 

T n = + . OI rw/1 I +0(/3 2 ,a 2 /3). (19) 

i?= (d 2 + 4d/3sm2fcL)/(l + d 2 ) J VM ' HJ K J 

Clearly unitarity, as conventionally defined, is violated in this process. 



3 Quasi-Hermitian Analysis 

The Hamiltonian has been specifically constructed to be PT symmetric. In 
such cases, as discussed in the introduction, we can in principle introduce 
a positive-definite metric operator[U [19] rj = e~® with respect to which H 
is quasi- Hermitian: W = rjHr]' 1 . Observables are those represented by 
quasi- Hermit ian operators A such that A^ = rjArj^ 1 . The original position 
operator x does not fall into this category: instead the observable of position 
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is X = pxp 1 , where p = 772 = e 2 Q . Consequently [201 I2l] . in this picture 
the relevant wave-function is not ip{x) = (x\i/j), but 

*(a?) = = <x|p|V> = y p{x,y)ip{y)dy. (20) 

Let us therefore attempt to construct $ for the present Hamiltonian. Un- 
fortunately this cannot be done exactly, but Mostafazadeh|18j has devised 
a general method for constructing a series expansion for rj in the coupling 
constants of a series of delta functions. This method is based on expressing 
the condition of quasi- Hermiticity of the Hamiltonian as a partial differential 
equation for rj(x,y), and converting it into an integral equation that can be 
solved iteratively. 

In some more detail, for the Hamiltonian H = p 2 + V(x), the integral 
equation for rj(x,y) takes the form 

rj(x, y) = u(x, y) + (Kr))(x, y), (21) 

where u(x, y) = u + (x — y) +u^(x + y) is the general solution of the differential 
equation (—9% + d^)u(x, y) — 0, and K is defined by 

/ drV{r) j ds+ dsV*(s) / dr J r?(s, r) (22) 

J J x—y+r J J y—x+s J 

In Ref. [18] the equation (T2TT) was written in the form 

ti = {1-K)- 1 u, (23) 

which in principle can be expanded in K, 1.6. clS du simultaneous expansion in 
the coefficients of the delta functions. However, in the present case we need to 
do something slightly different, because we are thinking of (3 as a perturbative 
parameter, but not a. Thus we need to split K up into K = K a + K^. Then 
we write Eq. (1211) as 

77 = u + (K a + Kp)r}. (24) 

Now for (3 = the Hamiltonian is Hermitian, so we want rj(x, y) = 5(x — y). 
This means that to order j3 we should take u = (1 — K a )8. In principle we 
could also add to u a term (3w, where w is a solution of the homogenous 
equation (d 2 — dy)w(x, y) = 0. However, it turns out that such a term is not 
required. So 

7} = (l-K a )6+(K a + K p )7i, (25) 
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so that 

(l-K a )r] = (1 - K a )5 + Kprj, (26) 

with solution 

77 = 6 + r^;^ 77 • (27) 

To O(P), which is as far as we will take the calculation, this reduces to 

1 = ^-rhrj^ 

= S + (l + K < , + K 2 c , + ...)K„6 . (28) 

Since we would like to treat (3 perturbatively, but not a, it is extremely fortu- 
nate that the higher powers of K a in this equation do not in fact contribute, 
as is shown in Section 3.2. 

3.1 Effect of Kp 

As Mostafazadeh has shown [TS] (with m = |, h — 1), the action of Kp on 
the delta function S(x — y) is 

(Kp6)(x, y) = ^i/3[9(x + y - 2L) - 6(x + y + 2L)\e{y - x), (29) 

where e(z) is the sign function e(z) = sgn(z). Thus, to order (3 we have 

[ x + y < -2L 

f](x,y) = 5(x — y) + < \i(3e{x — y) —2L<x + y<2L 

{ 2L < x + y 

= 6(x-y)-f3Q 1 (x,y) (30) 
3.1.1 Bound-state wave function 

Let us now calculate the effect of p(x,y) = 5(x — y) — \f3Qi(x,y) on the 
bound-state wave-function ty(x). From Eq. (1201) it is given, to this order, by 

i&(x) = ip(x) + -i(3 / e(x — y)ip(y)dy 



4 

i[)(x) - -i(3 



-2L-x 

2L~x rx 

ip{y)dy- / ip(y)dy 

J-2L-X 



(31) 



for < x < L. Here, since we are working to 0(0), we can take tp(y) = e K ^ 
in the integrands. Then 



= (1 + i(3e- 2KL )e- KX - i(3e- 2KL e KX 

»2L-x 



-HI] 



dy-2 



o 







/ e~ Ky dy- 


[ e Ky dy 


Jo 


'-2L-X 



(32) 



1 + ^i(3e~ 2KL \ e~ KX - U(3e~ 2KL e KX + 0(1 - e~ KX ) + 0(f3 2 ). 

The last term is not something we expect at all, but before jumping to 
conclusions we should await the calculation of the K a Kp contribution, which, 
because it involves a in the combination a = at/ (2k) = l+0(/3 2 ), is of exactly 
the same order. 

3.1.2 Scattering wave-function 

We are now concerned with the scattering wave-function for large \x\. For 
x > it is actually sufficient to take x > 3L (so that 2L — x< —L), at which 
point the wave-function, ^> settles down to its asymptotic form, namely 

1 /-2L-X 

= e ikx + -i(3 / e(x-y)[(l-ia)e iky + iae- iky ]dy 

4 J-2L-X 
i p2L-x 

= e ikx + -if3 [(1 -ia)e iky + iae~ iky ]dy (33) 

4 J-2L-X 

= (l-a(3 sin 2kL)e ikx + 0(1 - ia) sin 2kL e~ ikx + 0((3 2 ) . 
For x < —3L (so that — 2L — x > L) we instead get 



i rlL-x 

#<(a;) = Ae ikx + Be~ ikx + -if3 e(x-y)e iky dy 

4 J-2L-X 
1 p2L-x 

= Ae ikx + Be~ ikx - -i(3 / e iky dy (34) 

4 J-2L-X 

= (l-ia)e ikx + i[a + (3sm2kL~Aa(3sm 2 kL)]e- ikx + 0((3 2 ) . 

Note that the physical picture has completely changed, because we have 
incoming and outgoing waves on both sides. This is the rather drastic mod- 
ification noted in Ref. [16]. Some such modification at infinity is clearly 
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necessary if we are to restore unitarity in this picture. In a recent discretized 
model of scattering studied by Znojil[l7] the modification is instead a change 
in the normalization of the fluxes on either side of the scattering centre. In 
the present model we can check unitarity by comparing the net right-moving 
fluxes $ on each side. Unitarity is indeed restored to this order, because 

= = 1 - 2apsm2kL + 0(a 2 p,p 2 ). (35) 

Note that we are only allowed to calculate the fluxes using the standard 
formula in regions where the equivalent Hermitian Hamiltonian h is simply 
p 2 . Otherwise [1 6 j the conservation of probability takes a non-local form 
involving an integral over h(x,y). 

3.2 Effect of K a K p 

According to Eq. ([25]) we need to calculate K a Q x . In general[TS], the effect 
of K a on u(x, y) is 



(K a u)(x,y) = -~a 



■x+y 

9{y) I dru(r,0) + 9(x) / dsu(0,s) 

x—y J y—x 



(36) 



Recall, cf. Eq. ([30]), that Qi(x, y) = \ie{y -x)[9(x + y + 2L) - 9(x + y - 2L)\. 
Thus 

{K a Q x )(x,y) = -ia9{y) ! +V drs{r){9{r + 2L)-9{r-2L)} 

^ J x—y 

- jia9(x) ^ dse(s)[9(s + 2L) - 9{s -2L)]. (37) 

4 J y-x 

In principle we need to calculate the effects of K% arising from the expansion 
of Eq. (|28]) . However, a surprising and welcome result is that these vanish 
for n > 1. Thus, in the calculation of K 2 Qi according to Eq. (136]) . we need 



i r 

(KaQJfaO) = -jia9(r) / due(u)[9(u + 2L) - 9(u - 2L)} 

^ J — r 

= by symmetry. (38) 

Similarly (K a Qi)(0, s) = 0. So in fact all higher order terms in Eq. (128]) are 
absent. 
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3.2.1 Bound-state wave-function 



Here, since we are concerned with the limit as L — ► oo, we need to evaluate 
Eq. (l37j) for < i C L. A similar analysis will apply to i < 0, but this 
is easily obtained by PT symmetry. First it is easy to see that for a non- 
zero result x and y must have opposite signs. So in the present case y < 0. 
There are then three possibilities depending on the positions of the limits of 
integration y + x and y — x. The net result is 

I ( %y > y > -x 

(K a Qi)(x,y) = —ia I -2x -x > y > -2L + x (3,9) 

{ ~{2L + y + x) -2L + x > y > -2L - x 

giving a correction to the bound-state wave function 

A\l> = -ia(3 I 2 \ y -2x \ - (2L + y - x) \ e Ky dy 

° I J -x J-2L+X J-2L-X J 

= Ua(3[e- 2KL (e KX - e~ KX ) - 2(1 - e~ KX )] . (40) 

When added to ^ of Eq. fl32l) . and remembering that a — 1 + 0(j3 2 ), we 
obtain the remarkably simple result 

q(x) =e~ KX + 0{f3 2 ), (41) 

which to this order is equal to the original undisturbed wave function ijj(x) 
of Eq. ([2D for x > 0. 



3.2.2 Scattering wave-function 

For the corrections to the scattering wave-function we need to evaluate 
Eq. fl37j) for \x\ ^> L. Let us first consider i > L, in order to obtain 
the correction A^> to Again, for a non-zero result y must be negative, 
so that 

(K a Q 1 )(x,y) = --ia T ds e(s)[6(s + 2L) - 6(s - 2L)]. (42) 

^ J y—x 

Here the lower limit, y — x, is less than — 2L, and there are again two possi- 
bilities depending on the position of the upper limit. The net result is 

1 ( x + y-2L < x + y <2L 
(K a Q x )(x,y) = -ia < ~{x + y + 2L) -2L < x + y < , (43) 

y otherwise 
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giving 



--a/3 



dy(x + y - 2L)if> < (y) 



2L-x 
2L-x 



dy(x + y + 2L)il; < (y) 



where ij)<(y) — (1 — ia)e lky + iae %ky . The result of this calculation is 

A^> = 2d/3sin 2 kL [(1 - ia) e - ikx + iae ikx ] , 

giving the corrected value of as 

* > (a;) = (1 - a f3sin2kL + 2a 2 (3 sin 2 kL)e ikx 

+0(1 - id) (sin 2kL + 2a sin 2 kL)e~ ikx . 



(44) 



(45) 



(46) 



Now let us take x <C —L , so that y must be positive. The expression for 
K a Qi turns out to be the same as that given in Eq. f[59"j) . Thus 



A^< 



— -a/3 



dy(x + y - 2L)if>>(y) 



2L-x 
2L-x 



dy(x + y + 2L)il; > (y) 



where ip>(y) = e lky . Hence 

A^< = 2d/3sin 2 kLe~ ikx , 
giving the corrected value of \l/ < as 

*<(x) = (1 -ia)e ikx + i[a + (3 sin2kL-2af3 sin 2 kL)e~ ikx . 
From Eqs. (I4"6"j) and (I49p we obtain the fluxes 

$> = $< = !- 2a/3(sin 2kL - 2a sin 2 kL) + 0([3 2 ). 



(47) 



(48) 



(49) 



(50) 



Compared with Eq. (135|) . we now have the explicit 0(a 2 (3) terms, and the 
result is correct to the order shown. 
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4 Discussion 



Let us now consider the conceptual issues raised by these calculations. First 
it should be emphasized that we are concerned here with quasi-Hermitian 
Hamiltonians, that is, Hamiltonians that can be related by a similarity trans- 
formation to a Hermitian Hamiltonian. It is only for such Hamiltonians that 
we can attempt to construct the metric rj and to look at the situation in 
the quasi-Hermitian framework. For generic non-Hermitian Hamiltonians 
no such framework is available, and one is bound to treat them as effective 
Hamiltonians within the standard framework of quantum mechanics. In that 
case one would simply perform an analysis similar to that of Section 2 and 
accept that unitarity is not conserved, essentially because we are dealing with 
a subsystem of a larger system whose physics has not been taken fully into 
account. 

For the potential we have chosen we are able to consider both bound and 
scattering states. The bound state is the simpler to consider. The standard 
quantum mechanical analysis shows that the introduction of the perturbing 
PT-symmetric delta-function potentials does not significantly modify either 
the bound-state energy or the wave-function if these potentials are sufficiently 
distant. The analysis of Section 3, in particular the final result of Eq. (41), 
shows that this remains the case within the quasi-Hermitian framework, with 
the new bound-state wave-function ty(x) being identical to the original wave- 
function ip(x) for large L. This is a new and reassuring result, which we were 
unable to address in Ref. [16], where the potential did not support a bound 
state. The general message we would like to draw from this is that a localized 
physical Hermitian system is not significantly affected by the introduction 
of distant non-Hermitian potentials, and may be treated in the framework 
of standard quantum mechanics without the necessity of introducing a new 
metric. 

The real conceptual problems arise for the scattering states. The stan- 
dard quantum mechanical analysis shows that unitarity, as conventionally 
defined, is not conserved. In the quasi-Hermitian framework, however, one 
can hope that a modified form of unitarity is in fact conserved, as is indeed 
borne out by the calculations of Section 3. However, this involves a fairly 
drastic redefinition of asymptotic states, a non-local effect, given the finite 
support of the perturbing PT-symmetric potentials. The crux is the differ- 
ence between x, the coordinate parameter in terms of which the Hamiltonian 
is originally defined, and X, its quasi-Hermitian counterpart, defined by the 
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non-local relation X = pxp^ 1 . The former, x, is Hermitian and therefore 
an observable, in the standard framework of quantum mechanics, while X is 
not. Conversely, in the quasi- Hermitian picture, X is quasi-Hermitian and 
an observable, while x is not. The argument of the first wave-function ip(x) 
is the eigenvalue of x, while that of \fr(x) is the eigenvalue of X. 

What we have done in this paper is transform an initial scattering set- 
up defined in terms of x, and then consider the corresponding picture in 
the quasi-Hermitian framework in terms of X. The initial scattering set-up 
had plane waves entering from the left and then being either reflected or 
transmitted, with probability not being conserved. As we have seen, the 
corresponding quasi-Hermitian picture is that the newly-defined probability 
is indeed conserved, but that waves now enter from both left and right. An 
alternative mathematical possibility would be to set up a scattering situation 
in which ^> has only outgoing waves and then work backwards to construct 
t/>, which would undoubtedly have waves entering from both left and right. 
In either case, the physical picture changes drastically when going from one 
picture to the othefi 

In the author's opinion, the only satisfactory resolution of this dilemma 
is to treat the non-Hermitian scattering potential as an effective one, and 
work in the standard framework of quantum mechanics, accepting that this 
effective potential may well involve the loss of unitarity when attention is 
restricted to the quantum mechanical system itself and not its environment. 
This is indeed the attitude taken in a recent paper by Berry [23J , where var- 
ious intensity sum rules are derived for diffraction off PT-symmetric optical 
lattices. There it is taken for granted that the intensities are given by \ip\ 2 - It 
is true that these are classical calculations, but because of the correspondence 
principle the same thing would apply in quantum mechanics. 



Tin a recent paper, ref. [55], Znojil has constructed certain discrete matrix models of 
scattering in which in and out states are not mixed in the quasi-Hermitian picture, but 
instead the definition of the flux differs on the left and right of the scattering centre. This 
is a somewhat less drastic change, but still involves a departure from standard quantum 
mechanics at large distances. 
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